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Abstract. We extend the homogeneous model for the 
synchrotron emission from plerions to allow more realistic 
predictions of the surface brightness distribution at dif- 
ferent frequencies. Abandoning the assumption of a uni- 
form particle distribution, we assume that particles are 
injected into the synchrotron nebula only in the vicinity 
of the pulsar. Their distribution is then determined by 
MHD propagation in a magnetic field of spatially con- 
stant but time-dependent intensity. This highly simplified 
model reproduces the integrated spectrum and the syn- 
chrotron surface brightness profile of the Crab Nebula at 
radio frequencies. However, when applied to higher fre- 
quencies, it underestimates the extension of the emitting 
region, suggesting that particle diffusion with respect to 
the magnetic field lines must occur. 
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1. Introduction 

The integrated non-thermal spectrum of plcrionic super- 
nova remnants is well accounted for within the framework 
of the homogeneous model by Pacini & Salvati (1973, here- 
after PS), who regard the plerion as an adiabatically ex- 
panding spherical bubble of magnetized relativistic fluid. 
The bubble is replenished by the continuous conversion 
of the pulsar's energy outflow into magnetic energy and 
relativistic particles. The energy balance of the plerion 
is determined by the competition between the decreasing 
pulsar's input and the losses the fluid undergoes because 
of radiation and expansion. Both particles and magnetic 
field are assumed to be distributed uniformly through the 
bubble. 

When the synchrotron emission of the Crab Nebula is 
calculated on the basis of the PS model, both the decline 
of luminosity with time (Veron-Cetty & Woltjer 1991) and 
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the frequency of the low energy spectral break (Marsden 
et al. 1984) are consistent with a magnetic field strength 
around 3 x 1CT 4 G, in agreement with estimates based 
on equipartition arguments (Woltjer 1958) and also with 
recent measurements of the nebular inverse Compton radi- 
ation flux (De Jager & Harding 1992). Other applications 
of the PS model include the interpretation of the over- 
all spectral and evolutionary properties of other plcrionic 
supernova remnants, as 3C58 and Gil. 2-0. 3, which seem 
to be very different from the Crab Nebula (Woltjer et al. 
1997). 

Despite the success of the model at describing the in- 
tegrated properties of plerions, the assumption of homo- 
geneity is unsatisfactory. If as assumed the particles are 
accelerated in the general environment of the pulsar and 
if loops of magnetic field are created by its rotation, im- 
portant gradients could be expected in the particle and 
field distributions. Also observationally such gradients are 
in evidence. The Crab Nebula does not have a uniform 
emissivity at any wavelength: it is concentrated towards 
the center, the more strongly so at the higher frequencies. 

We shall here consider a modification of the PS model 
in which loops of magnetic field are continuously injected 
with the relativistic electrons attached to the field lines. 
Since little is known about the detailed conditions around 
the pulsar we shall assume that the injection takes place 
at a certain distance r$ from the center. In the case of the 
Crab Nebula the moving wisps do in fact suggest that the 
injection occurs at some 10" from the pulsar. Since the 
propagation velocity of magnetic disturbancies is likely to 
be much larger than the expansion velocity of the bubble, 
we shall assume that the field adjusts itself to constant 
energy density, an assumption that is strictly justified only 
if the particle pressure is below the magnetic pressure and 
the field is chaotic. 

While in the case of the Crab Nebula the polarization 
studies at radio and optical wavelengths (Wilson 1972a, 
Woltjer 1958) indicate a predominance of azimuthal mag- 
netic field, the modest degree of polarization shows that 
there are different field directions along the line of sight. 
Such a less regular field (and the effect of reconnection) 
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may well cause the particles to move more freely through 
parts of the nebula. Thus the present model, which in- 
cludes only convection, is likely to be a limiting case. In 
fact we will find that the observed emissivity distribution 
is intermediate between the one computed here and the 
homogeneous case. 

2. Outline of the model 

In our model the nebula is regarded as a sphere of radius 
R(t), expanding at a constant rate. A continuous supply 
of new magnetic energy and particles is provided by the 
central pulsar, whose energy output per unit time can be 
written as 



L(t) 



(l + i/VK 



(1) 



Both new particles and new magnetic flux are injected 
into the nebula within a spherical layer at a distance r 
from the central pulsar. We shall assume that the injection 
region has zero extension, although very likely this is an 
oversimplification. 

We assume that the particles' velocity distribution is 
isotropic and that each particle emits essentially at its 
characteristic frequency: v c — c 2 B±E 2 , where c 2 = 0.29x 
(3ec/47r(mc 2 ) 3 ) is a constant, B± = Bain9, and E and 9 
are the particle's energy and pitch angle, respectively. Un- 
der these assumptions the nebular synchrotron emissivity 
becomes 
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with Ci = (2c/3)(e/mc 2 ) 4 a constant and N(E, t, r) repre- 
senting the particles' spectral and spatial density at time 
t. The integral represents the average over particles' pitch 
angles. 

As stated before, we regard the magnetic field as spa- 
tially constant, therefore, in our model, both the spatial 
and spectral characteristics of the nebular synchrotron 
emission are directly related to the particle spatial and 
spectral distribution, N(E,t,r). 

N(E 7 t, r) is related to the injected particles' spectrum 
through the particles' number conservation law: 



N(E,t,r) 
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with J being the number of newly injected particles per 
unit time and energy interval and Ei and U being, re- 
spectively, the initial energy and injection time of a given 
particle. The last term in Eq. (3) represents the jacobian 
of the transformation (ti,Ei) — > (r, E). 

Apart from knowledge of the injected spectrum, in or- 
der to calculate N(E,t,r) from this equation, knowledge 
of the particles' age (t — U) and initial energy (Ei) as 
functions of their present position (r) and energy (E) is 
needed. 



We consider the energy evolution of each particle due 
to adiabatic and synchrotron losses: 



dE_ 
~dt 



E. 

3 



= -=-V -u- Cl B 2 ± E 2 . 



(4) 



Once the magnetic field B and the fluid bulk velocity u are 
known at each time, Eq. (4) can be integrated to determine 
Ei(E,r). 

The determination of the time dependence of the mag- 
netic field strength B is straightforward under our as- 
sumptions of time constancy of the nebular expansion rate 
and spatial constancy of B itself. Following PS, we re- 
late B to the nebular content of magnetic energy Ws(t): 
B 2 (t)=8 ir W B (t)/ V(t), with V(t) = (4ir/3)(R 3 (t) - r 3 ) 
representing the confining volume for the magnetic field. 
Then we write the time evolution of Wb due to expansion 
losses and the decreasing pulsar input: 



dW B W B d 

~ = " — ^ lnT/(<)+/3i(<) • 



(5) 



Integration of Eq. (5), where [3 is a constant, yields the 
magnetic field strength at each time (see, for instance, PS) 

Knowing B(t), and further assuming that the magnetic 
field is largely azimuthal, B = B(t)e<p, we can use the flux 
freezing condition, 



dB „ , 

— = VA (uAB) 

dt v ' 
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as an equation for the fluid bulk velocity u, which we 
assume to be radial: u = u(r,t)e r . Integrating Eq. (6) 
with the boundary condition that the fluid velocity field 
matches the expansion velocity of the nebula at its outer 
edge (u(R(t),t) —v), we find: 



u(r,t) = 



vR(t) 8\nB (R 2 (t)-r 2 ) 



+ 



dt 



2r 



(7) 



The velocity field starts at r at a fraction of the speed of 
light of order (v/c) x (R/ro) (~ 1/10, see below) and then 
it decreases roughly as 1 jr to match the nebular expansion 
velocity at R. 

Eq. (7) can also be written as 



dt 



[B(t)(R 2 (t)-r 2 (t))}=0, 



(8) 



which yields, in its integral form: 



B(t)(r 2 -r a 2 )= (9) 

B(t) (R 2 (t) - r 2 ) - B(U) (R 2 (U) - r 2 ) . 

This last equation simply states that the magnetic flux 
injected into the nebula during the time interval t — tj, 
with U being the time at which a particle that at time t is 
at r was born, is all contained between r and r. Taking 
r = r(t,ti) from Eq. (9), we have got the relation needed 
to connect the particles' position at each time to their age. 
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Presently what is left to determine, before being able to 
calculate N(E, t, r) as a function of the injected spectrum, 
is the evolution of the particles' energy from initial to 
present value. Introducing the expression found for u into 
Eq. (4), we obtain: 



d 
dt 
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We integrate Eq. (10) after replacing sin 2 # with its time 
average. We take this to be 2/3, which is equivalent to 
state that each particle during its synchrotron lifetime ex- 
periences all the possible velocity orientations with respect 
to the magnetic field with the same probability. Obviously, 
we expect this approximation to work better the longer a 
particle lives, hence the smaller its initial energy is. Nev- 
ertheless we apply it to particles of all energies and finally 
find the expression for Ei(E, t, r): 
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The meaning of the energy Eb(t,r) is apparent: it repre- 
sents the maximum possible energy for particles that at 
time t reside at r. Substituting £& with yj t//(c2-Bj_), Eq. 
(11) also defines the maximum radial distance rM(v,t) 
from which we expect emission to come for a given fre- 
quency. 

Inserting our findings into Eq. (3) we are finally able 
to relate N(E, t, r) to the injected particle spectrum. Con- 
cerning the latter some assumptions are necessary. As- 
suming that the injected electron spectrum is described 
by a single power law, synchrotron aging can account for 
just one break in the emission spectrum, while observa- 
tion shows that in the case of the Crab Nebula there is 
another break between the optical and the X-rays; this 
break must be intrinsic to the injection mechanism. Then 
we allow for a particle energy spectrum with two different 
slopes in two different energy ranges, namely: 



K^t) E-^ ■ E m < E < Ei 
J{E,t) = { K 2 {t) E-^ ; E X <E< E M 
; E < E m or E > E M ■ 



(13) 



For simplicity, we assume all the energy cuts to be 
time independent and the ratio between K\ and K 2 to 
remain constant. Since we have assumed that a constant 
fraction [3 of the pulsar's energy outflow goes into feeding 
the magnetic field, the remaining 1 — (3 will be used for 
accelerating particles, and then the constraint on the total 
slowing-down power converted into particles, 



(1 



-(3)L(t)= / 



J(E,t)EdE, 



(14) 



implies 

Kt(t) cx K 2 (t) oc L(t). 



(15) 



Finally we have for the spectral and spatial distribu- 
tion of particles across the nebula: 
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71, E m (r) < E < E^r) 

72, £i(r) < E < E M {r) 



(16) 



(17) 



where the energies E m (r), Em{t) and E\{r) are the 
evolved minimum, maximum and intrinsic break energy 
at radius r, respectively, and can be calculated as func- 
tions of E m , Em and E\ by means of Eq. (11). 

3. The case of the Crab Nebula 

For the Crab Nebula pulsar we have in Eq. (1) Lo = 
3 x 10 39 erg/s, r = 710 yr, n = (n + l)/(n - 1) where 
the braking index n is n — 2.5 (Groth 1975). For the pa- 
rameters of the injected spectrum (Eq. (13)), following a 
procedure similar to PS, we derive from the observations 
the values 71 = 1.54, 72 = 2.3, E 1 ~ 20eV. The parame- 
ter P in Eq. (14) determines, through Eq. (5), the present 
value of B; we take the latter to be equal to 3 x 10~ 4 G, so 
that — 0.25. Given these parameters the only unknowns 
in Eq. (16) are the injection radius r and N, the latter 
containing the cut-off energies of the injected spectrum. 
These parameters have very different observational signa- 
tures: the first one is related to the radial distance from 
the central pulsar of the luminosity peak, and the second 
one only affects the overall nebular synchrotron flux. 

The model predictions for the synchrotron surface 
brightness profile of the Crab Nebula have been compared 
with high resolution data at various frequencies. 

For the radio band we have used a VLA map at a 
frequency of 1.4 GHz (Bietenholz et al. 1997). 

The spatial and spectral distribution of the optical 
synchrotron continuum was determined by Veron-Cetty 
& Woltjer (1993) after subtraction of the thermal con- 
tributions from foreground stars and filaments, from four 
narrow-band images at wavelcnghts of 9241, 6450, 5364 
and 3808 A. We have reanalysed these maps, kindly put 
at our disposal by M.P. Veron-Cetty, with state-of-the-art 
star subtraction algorithms. 

Finally, in the X-ray band, we have used, after de- 
convolution of the instrumental PSF and subtraction of 
the dust halo (Bandiera et al. 1998), a collection of all the 
public ROSAT HRI data concerning the Crab Nebula. We 
have estimated the mean photon energy of these data to 
be 1 keV. 
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In our model the synchrotron surface brightness J„ is 
simply obtained by integration of Eq. (2) along the line of 
sight: 

f m rdr 
Mv,t,z)= S„{is,t,r)-==, (18) 

J z \ v z 

with the expression for the particle number density N in 
S v calculated from Eq. (16). 

In order to compare our spherical model with the ob- 
servations, we have extracted from each image what we 
call a "radial intensity profile" : we first sampled the emis- 
sion profiles of the nebula along different directions, taking 
the mean values over small areas of 10" x 10" and then, 
after rescaling the different axes to a common length, we 
averaged those profiles. The procedure just described, to 
which we refer in the following as "data sphericization" , 
is the main cause of uncertainty in our radial profiles and 
it is what we take into account in the error bars attached 
to the data points in most of the following plots. 

When the radiative losses arc negligible (E/E b (r, t) <C 
1 in Eq. (16) throughout the entire nebula) N and B 
enter the expression of the synchrotron emissivity sim- 
ply as multiplying factors. Therefore fitting the shape of 
the surface brightness profile at radio frequencies allows 
a straightforward determination of ro, and once this is 
known, the value of N simply comes from fitting the inte- 
grated flux, so that the model is fully determined. 

Our best fit estimate gave ro = 10", which, for a 2 
kpc distance to the Crab Nebula, translates into ro = 0.1 
pc. The corresponding radial profile at 1.4 GHz is plotted 
against the data in Fig. 1. This value of r yields a distance 
from the pulsar to the injection site fully compatible with 
the association between the particle acceleration region 
and the location of the optical wisps. 
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Fig. 1. Comparison with the radio data (points) of the PS 
model (dotted curve) and of our best fit model (solid curve) 
for the surface brightness profile of the Crab Nebula. The ab- 
scissa is the projected distance from the pulsar. The error bars 
attached to the data points take into account the uncertainties 
introduced by "data sphericization" . 



As in the PS model, the integrated fluxes reproduce 
the observations at all frequencies: both the solid and the 
dashed curve yield the same flux as the interpolation of 
the data, when integrated on a spherical surface of ra- 
dius 2 pc. Nevertheless, although in the radio part of the 
spectrum the fit to the observed profile is rather good and 
substantially improves the homogeneous model, the model 
predictions fail to reproduce the data at optical and X-ray 
wavelengths. 
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Fig. 2. Comparison with the data at radio (1.4 GHz), opti- 
cal (6450 A) and X-ray (1 keV photons) frequencies of the 
expected size of the emitting region. We plot as a function 
of frequency the radius within which half of the total nebular 
flux per unit frequency is produced. The solid curve is based 
on the present model whereas the dotted curve represents the 
homogeneous PS model. Again the error bars attached to the 
data points take into account the effect of our "sphericization" 
procedure. 

As shown in Fig. 2 the emission at optical and X-ray 
wavelengths calculated on the basis of the present model is 
too concentrated: the highest energy particles emit most of 
their energy immediately after the injection. This causes 
most of the flux to originate from a narrow region and the 
particles to travel a very short distance from the injection 
site before their energy is degraded by severe synchrotron 
losses. 

This effect could be cured by substantially lowering the 
magnetic field and thereby the synchrotron losses. How- 
ever the inverse Compton data do not allow this. More- 
over if the magnetic field energy is less than the particle 
energy the model becomes invalid and at lower fields the 
total particle energy would soon exceed that produced by 
the pulsar. 

The homogeneous PS model in which the particles 
move freely through the Nebula yields too broad a dis- 
tribution of the emissivity, our model a too narrow one at 
the higher frequencies. Apparently some of the particles 
can reach the outer parts of the nebula without suffering 
the large synchrotron losses which occur when they are 
fully tied to the field lines. 
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Diffusion of particles could solve this problem, but the 
diffusion coefficient would have to be ~ 10 times larger 
than for Bohm diffusion, in agreement with previous esti- 
mates (Wilson 1972b). A more complex field structure in 
which some lines connect the inner and outer parts of the 
nebula might give an acceptable solution, since the parti- 
cles could move along these lines at a substantial fraction 
of the speed of light. 
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